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ABSTRACT 

We employ high-spatial resolution spectropolarimetric observations from the Solar Optical Telescope on-board the Hinode spacecraft 
to investigate the fine structure of the penumbral magnetic fields. The Stokes vector of two neutral iron lines at 630 nm is inverted at 
every spatial pixel to retrieve the depth-dependence of the magnetic field vector, line-of-sight velocity and thermodynamic parameters. 
We show that the azimuthal angle of the magnetic field vector has opposite sign on both sides above the penumbral filaments. This 
is consistent with the wrapping of an inclined field around the horizontal filaments. The wrapping effect is stronger for filaments 
with larger horizontal extensions. In addition, we find that the external magnetic field can penetrate into the intraspines, leading 
to non-radial magnetic fields inside them. These findings shed some light on the controversial small-scale structure of the sunspot 
penumbra. 
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1. Introduction 

Our knowledge of the small-scale structure of magnetic features 
in the solar photosphere remains hindered by the limited spatial 
resolution we can achieve. A prominent example of this situation 
is the fine structure of the sunspot penumbra, where a variety of 
different models compete to explain the observations (Solanki & 
Montavon 1993; Schlichenmaier et al. 1998; Sanchez Almeida 
2005; Spruit & Scharmer 2006; Borrero 2007). The degree of 
sophistication of these models is such that we can no longer use 
moderate spatial resolution (~ 1") observations to distinguish 
among them. 

Many of these models assume the presence of two different 
magnetic components: one in the form of a somewhat inclined 
(~ 40 - 50° with respect to the vertical) and strong (~ 2000 G) 
magnetic field, and another in the form of a weaker and more 
horizontal magnetic field. They are usually referred to as spines 
and intraspines, respectively. Their presence has been repeat- 
edly confirmed observationally (Lites et al. 1993; Riiedi et al. 
1998; Borrero et al. 2004,2005; Bellot Rubio et al. 2004; Bello 
Gonzalez et al. 2005; Langhans et al. 2005). 

Some of these models identify the intraspines (horizontal 
and weak magnetic field) with a horizontal flux tube (uncombed 
model), whereas in others it appears as the consequence of a 
field free gap protruding from beneath [gappy model). Either 
way, spines are always assumed to avoid the intraspines by 
folding around them (Martinez Pillet 2000). Although stud- 
ies of the azimuthal variation of the Net Circular Polarization 
(Schlichenmaier et al. 2002; Miiller et al. 2002; Borrero et al. 
2007) and spectropolarimetric observations at moderate spatial 
resolution (Lites et al. 1993) provide an indirect measure for this 
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effect, direct observational evidence supporting this assumption 
has never been presented. 

In this work we use Hinode's high-spectral [AA/A < 4xl0~ 6 ) 
and high-spatial (~ 0.32") resolution polarimetric observations 
of the pair of photospheric Fe I lines at 630 nm to retrieve the 
depth dependence of the magnetic field vector. In section 2 we 
describe those observations. Section 3 presents the inversion pro- 
cedure that allows us to infer the magnetic field vector from the 
recorded polarized spectra. Section 4 presents our results and fi- 
nally, Section 5, the conclusions. 

2. Observations 

On May 3rd 2007, between 10:15 and 11:40 am UT, the active 
region AR 10953 was mapped using the spectropolarimeter of 
the Solar Optical Telescope on-board of the Hinode spacecraft 
(Lites et al. 2001, Kosugi et al. 2007; Shimizu et al. 2007). The 
active region was located at a heliocentric angle of 8 = 19.2°. It 
was scanned in a thousand steps, with a step width of 0. 148" and 
a slit's width of 0.158". The spectropolarimeter recorded the full 
Stokes vector (/, Q, U and V) of the pair of neutral iron lines at 
630 nm with a spectral sampling of 21.53 mA. The integration 
time was 4.8 seconds, resulting on an approximate noise level of 
1.2 x 10 3 . In the absence of the telluric oxygen lines we pro- 
ceeded with two different wavelength calibration methods that 
were cross-checked for consistency. The first method was ob- 
tained by matching the average quiet Sun profile with the FTS 
spectrum, whereas the second calibration assumes that the av- 
erage umbral profile exhibits no velocities, which according to 
Rezaei et al. (2006) is a reasonable assumption. 

Map of the continuum intensity on the limb-side at 630 nm 
is shown in Figure 1 . The black arrow indicates the direction of 
the center of the solar disk. The penumbra on the center side is 
heavily distorted and therefore left out from our analysis. On the 
limb side the penumbra is more uniform, with radially aligned 



2 



Borrero et al,: Magnetic field wrapping in penumbral filaments 



AR 10953 yiA=0.866 



30 

25 

20 

15 

10 

5 








' I I I I I I I I I 


I Y : 

1 1 1 1 1 1 1 1 1 1 I 1 1 1 1 I 1 1 1 1 1 



10 15 20 
X [arcsec] 



25 30 



Fig. 1. Continuum intensity map at 630 nm of the limb-side of 
AR 10953. The umbra has been removed to increase the image 
contrast. This sunspot was observed using Hinode's spectropo- 
larimeter on the 3rd of May, 2007 at an heliocentric angle of 
= 19°. The black arrow points towards the center of the solar 
disk. The pair of dashed lines indicated by Y* lie near the line 
of symmetry of the sunpot. Results along these two slices are 
discussed in Section 4. 



filaments. This is the region that we have chosen to study. This 
sunspot has negative polarity (magnetic field in the umbra points 
towards the solar interior), however the results presented here- 
after are shown, in order to facilitate the interpretation, as if the 
sunspot had positive polarity. 

3. Inversion procedure 

To retrieve the physical parameters of the solar atmosphere from 
the spectropolarimetric observations, we employ the SIR code 
(Ruiz Cobo & del Toro Iniesta 1992). This code allows all rel- 
evant physical parameters to be a generic function of the opti- 
cal depth: T(f), B{t),j{t), <p(r), Vi os (t), etc. SIR retrieves the 
values of those parameters at a number of optical depth points 
called nodes. The final stratification is obtained by interpolat- 
ing splines across those nodes. Given the high-spatial resolution 
of our observations we assume that the penumbral structure is 
horizontally resolved, thus we perform an inversion with only 
one magnetic component, for which we allow 5 nodes in T(t), 2 
for B(t), y(r), (p(r) and Vi s(t), 1 for V mac and V m ; c (macro and 
microturbulent velocities). Including an extra free parameter to 
account for the scattered light, we have a total of 16 free pa- 
rameters. The same scattered light profile is used in all inverted 
penumbral pixels. It is obtained by averaging the intensity pro- 
files of those pixels with polarization signals below the noise 
level. The inversion retrieves typical values for the amount of 
scattered light of 8 to 17 %. 

As a result of applying the SIR inversion code, we obtain the 
magnetic field vector in the observer's reference frame at each 
spatial pixel as a function of the optical depth. Q To facilitate the 



1 An inversion of Hinode data from a different sunspot has been pre- 
sented by Bellot Rubio et al. 2007 and Ichimoto et al. 2007. They con- 
sidered only the Milne-Eddington case in which physical parameters 
are assumed to be height-independent. 



interpretation of the results it is convenient to transform to the 
local reference frame. In this paper we will discuss only those 
results obtained along the inclined dashed lines in Fig. 1 (de- 
noted as Y*). This region lies close to the line-of-symmetry of 
the sunspot (the line that connects the center of the sunspot with 
the center of the solar disk). In this region the conversion to the 
local reference frame, X*Y*Z*, can be done as follows: 



B*Jt) 



B(r)[siny(r) cos 0(t) cos 9 - cosy(r) sin 6] 
B{t) siny(r) sin 4>(j) 
B*(t) = B(t) [cos y(r) cos 6 + sin y(r) cos 0(t) sin 6] 



b;,(t) 



(i) 

(2) 
(3) 



where 8 corresponds to the heliocentric angle of the sunspot and 
B(t), j{t) and 0(t) are obtained from the inversion. Note that the 
knowledge of 0(t) is affected by the inherent 180 degrees ambi- 
guity. In order to distinguish between the two possible solutions: 
<p and <p + n, we consider that the magnetic field must point out- 
wards from the sunspot center and therefore B* x > 0. This allows 
one to uniquely determine the azimuthal angle. Finally, the in- 
clination and azimuth of the magnetic field vector in the local 
reference frame may be obtained through: 



*F = tan 1 

£ = cos 



(4) 
(5) 



4. Results and discussion 

Figures 2 and 3 shows the vertical distribution of the line-of- 
sight velocity (upper left panel), magnetic field strength (upper 
right panel), magnetic field inclination £ (lower left panel) and 
magnetic field azimuth *F (lower right panel), along the left and 
right slices Y* respectively (dashed lines in Fig. 1). Spinal and 
intraspinal structures are clearly visible. These examples dis- 
play 6 full intraspines, denoted as il through i6 in the velocity 
maps. Their vertical extension is about 1-1.5 logarithmic units 
of optical depth (roughly 150-250 Km). The horizontal exten- 
sion shows larger variations, ranging from 300 (i2) to 1000 Km 
(il). This is consistent with the large pixel-to-pixel variations in 
the filling factor between the strong/vertical and weak/horizontal 
components obtained by Bellot Rubio et al. (2004), using a two- 
component inversion of spectropolarimetric observations at 1" 
resolution. From the inversion of other regions in the penum- 
bra (not shown here) we also see examples of some very large 
intraspines that appear to be formed by two smaller ones very 
close to each other. 

The magnetic field shows a distinct pattern in azimuthal an- 
gle. Above the intraspines ¥ changes sign, being negative on the 
upper-left region above the intraspines, but positive on its right. 
This effect is highlighted by the thick solid line in the maps of 
the line-of-sight velocity, field inclination and azimuth, which 
show the variations of T* at an optical depth level of t = 0.1. 
This effect unmistakably denotes a change in the sign of B* and 
therefore indicates that the spinal magnetic field wraps around 
the intraspine. This is further confirmed by over-plotting the vec- 
tor field (see arrows in Fig 2 and 3) in the Y*t plane0. Note that 



2 We do not convert to vertical scale Z* from the optical depth scale r 
since this requires additional assumptions (i.e.: hydrostatic equilibrium) 
that are not fully justified. 
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B* x is missing in these figures, therefore the real magnetic field 
vectors will be pointing outside of that plane. 

The wrapping effect is seen in four of the six presented ex- 
amples. Above the two smaller intraspines (i2 and i5) *F is neg- 
ative on both sides, although it comes close to being positive 
on the right. This might indicate that horizontally narrow in- 
traspines perturb the surrounding field somewhat less than thick 
intraspines, where the external field is strongly forced to bend 
around them. 

Lites et al. (1993) also found that ¥ has opposite signs on 
both sides of the intraspines. They carried out a Milne-Eddington 
inversion of ASP data at lower spatial resolution. In our study we 
used a more sophisticated inversion technique that accounts for 
the optical depth dependence of the physical parameters. This 
kind of inversion is particularly appropriate in the limb side of 
the penumbra, where highly asymmetric Stokes profiles (often- 
times multi-lobed Stokes V) are clear indicative of gradients 
along the line-of-sight in the magnetic field and velocity. As a 
result, we have also detected that the magnetic field in the spines 
sometimes avoids the intraspines. This is the case of i2 and i6, 
where at optical-depth levels t 5 e [1, 0.1], the magnetic field has 
very small or no B* component: *f - and therefore B B x * 
(magnetic field inside the intraspines is aligned with the radial 
direction in the penumbra). In the rest of the examples: il, i3, i4 
and i5, due to the strong non-vanishing B* component, the mag- 
netic field inside the intraspines does not point radially outwards 
from the sunspot. Borrero (2007) has presented a magnetohydro- 
static model for cylindrical flux tubes where the magnetic field 
is not aligned with the tube's axis. Although our results cannot 
be used to confirm or rule out that model, an analysis targeting 
this particular effect, could be more conclusive. 

The Evershed flow is mainly concentrated in the intraspines, 
with velocities reaching up to 4.5 km s _1 . Velocities towards the 
observer (blueshifted velocities; Vi os < 0) appear sometimes 
high above the intraspines. This feature could bear important 
similarities to the upflows seen above umbral dots in 3D MHD 
simulations (Schiissler & Vogler 2006), however that possibil- 
ity is ruled out by the fact that inversions of our observations 
in the center-side reveal red-shifted velocities in the top layers 
of the intraspines (which are now characterized by blue-shifts 
in their lower layers). Therefore this pattern seems more com- 
patible with the Evershed & inverse Evershed flow (Deming et 
al. 1988; see also Fig. 11 in Bellot Rubio et al. 2006). Another 
possibility is that this effect is only an artifact of the inversion 
process. In deep layers we observe strong red-shifted velocities, 
whose magnitude decays rapidly as we move towards higher lay- 
ers. Since we only allow for 2 nodes in velocity (linear behavior 
with optical depth) the resulting velocities at t 5 ~ 10~ 3 could 
be blue-shifted simply because no other option is allowed to the 
inversion code. We have repeated our inversions using 3 and 4 
nodes and upflows, although weaker in magnitude, could still be 
seen above the intraspines. In the light of these results we cannot 
draw a decisive conclusion, but it is certainly something worth 
studying in the future. 

5. Conclusions 

Using high spatial resolution spectropolarimetric observations 
recorded by the Hinode spacecraft, we have demonstrated for 
the first time that the magnetic field in the penumbral spines 
(strong and vertical magnetic field) folds and bends around the 
intraspines (weaker and more horizontal magnetic field) while 
avoiding them. It is also frequent to observe how the trans- 
verse component of the external magnetic field leaks into the in- 



traspines. The overall magnetic configuration is at odds with the 
MISMA model (Sanchez Almeida 2005) for the penumbral mag- 
netic field. However, it partially supports the geometry adopted 
by the uncombed model (Solanki & Montavon 1993, cf. Borrero 
2007) and the gappy model (Spruit & Scharmer 2006). 

Although at this point our results seem to match better with 
the uncombed penumbra, the agreement is not perfect. In addi- 
tion, the highly simplified version of the gappy penumbra pre- 
sented in Scharmer & Spruit (2006) could be modified into a 
configuration similar to that of Schiissler & Vogler (2006), that 
presents more similarity to our findings. As our next step we will 
carry a more thorough analysis of the differences and agreements 
between models and observations. 

Acknowledgements. Thanks to the referee, Rolf Schlichenmaier, for his useful 
comments. In particular one that lead to discover an error in our calculations. 
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Fig. 2. Two dimensional cut across the leftmost Y* slice in Fig. 1 showing the depth-dependence of the line of sight velocity (upper- 
left panel), magnetic field strength (upper-right panel), magnetic field inclination in the local reference frame ( (lower-left panel) 
and magnetic field azimuth in the local reference frame ¥ (lower-right panel). The magnetic field vector at different locations in the 
Y*t plane is indicated by the arrows. 
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Fig. 3. Same as in Figure 2 but for the rightmost slice across Y* in Figure 1. 



